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ABSTRACT 

We investigate environmental effects on evolution of bright cluster galaxies (L > L,) in a A dominated 
cold dark matter universe using a combination of dissipationless A-body simulations and a semi-analytic 
galaxy formation model. The A-body simulations enable us to calculate orbits of galaxies in simulated 
clusters. Therefore we can incorporate stripping of cold gas from galactic disks by ram pressure (RP) 
from intracluster medium into our model. In this paper we study how ram pressure stripping (RPS) and 
small starburst induced by a minor merger affect colors, star formation rates (SFRs), and morphologies of 
cluster galaxies. These processes are new ingredients in our model and have not been studied sufficiently. 
We find that the RPS is not important for colors and SFRs of galaxies in the cluster core if star formation 
time-scale is properly chosen, because the star formation is sufficiently suppressed by consumption of 
the cold gas in the disks. Then observed color and SFR gradients can be reproduced without the RPS. 
The small starburst triggered by a minor merger hardly affects the SFRs and colors of the galaxies as 
well. We also examine whether these two processes can resolve the known problem that the hierarchical 
clustering models based on the major merger-driven bulge formation scenario predict too few galaxies 
of intermediate bulge-to-total luminosity ratio [B /T) in clusters. When the minor burst is taken into 
account, the intermediate B /T population is increased and the observed morphology gradients in clusters 
are successfully reproduced. Without the minor burst, the RPS cannot increase the intermediate B/T 
population. On the other hand, When the minor burst is considered, the RPS also plays an important 
role in formation of the intermediate B/T galaxies. We present redshift evolution of morphological 
fractions predicted by our models. The predicted number ratios of the intermediate B/T galaxies to 
the bulge-dominated galaxies show nearly ffat or slightly increasing trends with increasing redshift. We 
conclude that these trends are inevitable when bulges are formed through mergers. We discuss whether 
our results conflict with observationally suggested Aso/Ae evolution in clusters, which is a decreasing 
function of redshift. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: clusters: general — galaxies: 
halos — galaxies: interactions 



1. INTRODUCTION 

Rich clusters of galaxies are the most suitable labo- 
ratories to study how galaxy environments affect galaxy 
properties. The pioneering work by Butcher & Oemler 
(1978, 1984) has demonstrated that distant clusters con- 
tain much higher fraction of blue galaxies than nearby clus- 
ters, though it has been questioned by Andreon & Ettori 
(1999). It has been also found that galaxy morphology is 
a function of its environment (Dressier 1980; Whitmore, 
Gilmore, & Jones 1993). The most recent advances have 
been made with the Hubble Space Telescope, which allows 
the morphology of distant galaxies to be directly compared 
with the properties of their nearby counterparts. Dressier 
et al. (1997) and Couch et al. (1998) have suggested that 
the predominant evolutionary effects are that the distant 
clusters have a substantial deficit of SO systems compared 
with nearby clusters, and at lower luminosities they con- 
tain substantial Sc-Sd spirals, compared with the large 
population of dwarf spheroidals in present-day clusters, 
while it has also been questioned by Andreon (1998). On 



the other hand, some recent studies (e.g. Abraham et al. 
1996; van Dokkum et al. 1998; Balogh et al. 1999) have 
focused on observed trends in star formation rate (SFR) 
and morphology as a function of position within a cluster. 
These studies have shown that there is smooth transition 
from a blue, disk-dominated population of galaxies in the 
outskirts of clusters to a red, bulge-dominated population 
in the cluster cores. 

Many authors have suggested that the predominance 
of red, early-type galaxies in local clusters is the result 
of mechanisms that suppress star formation in high den- 
sity environments. This suppression makes galaxies within 
cluster cores redder and leads to a transformation of galaxy 
morphology. Comparison between the galaxy populations 
of local and distant clusters or those in cluster cores and 
outskirts provides a strong evidence for this scenario. This 
leads to a naive conclusion that the primary effect of the 
cluster environment is to transform luminous spiral galax- 
ies into SO galaxies through the suppression of their star 
formation. 
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Several mechanisms that may suppress the star forma- 
tion and transform one morphological type into another 
have been proposed. 

Interactions between galaxies are one possible process to 
promote morphological transformation. A'^-body simula- 
tions confirmed that major mergers between disk galaxies 
produce galaxies resembling ellipticals as merger remnants 
(e.g. Barnes 1996) and that accretion of small satellites 
onto their host spiral lead the morphology of the host spi- 
ral to SO type (Walker, Mihos, & Hernquist 1996). The 
galaxy mergers trigger the starburst, therefore the cold 
gas contained in these galaxies are exhausted in a very 
short time. These processes then cause the truncation of 
star formation. Although the bulge formation scenario 
by major mergers gives an excellent explanation both for 
colors (Kauffmann & Chariot 1998; Nagashima & Gouda 
2001) and distribution of cluster ellipticals (Okamoto & 
Nagashima 2001; Diafcrio et al. 2001), the observed SO 
population cannot be reproduced only by major mergers 
(Okamoto & Nagashima 2001; Diaferio et al. 2001). In ad- 
dition, the galaxy mergers are predominant effect before 
massive cluster formation, because, in clusters, the rela- 
tive velocity of galaxies is too high for such mergers to 
be frequent (Ghigna et al. 1998; Okamoto & Habe 1999). 
Therefore the galaxy merger is not promising as a key 
mechanism behind the Butcher-Oemler effect. Moore et al. 
(1996) examined the effects of rapid gravitational encoun- 
ters between galaxies. This mechanism has been called 
galaxy harassment and is highly effective at transforming 
fainter Sc-Sd galaxies to dSphs. While the galaxy harass- 
ment can account for the observed evolution of faint galax- 
ies in clusters, the concentrated potentials of luminous 
Sa-Sb galaxies help to maintain their stability, then the 
bulge-to-disk luminosity ratios {B/Ts) of these galaxies 
are hardly changed only by the galaxy harassment (Moore 
et al. 1999). Nevertheless, if the truncation of star for- 
mation in these galaxies takes place, these galaxies can 
become the objects that are similar to SO's, because the 
disks of Sa-Sb galaxies are substantially thickened and the 
spiral features are vanished by the harassment. It should 
be noted that correlation between B/T and Hubble type 
has significant scatter (Baugh, Cole, & Frenk 1996) and 
there is no guarantee that this relation is applicable to 
cluster galaxies. Accordingly we adopt the classification 
based on B/T to observed galaxies as well as to our model 
galaxies. 

The second mechanism is the truncation of star forma- 
tion through the removal of the diffuse hot gas reservoir 
that is confined in galactic dark halos and surrounds galax- 
ies (Larson, Tinsley, & Caldwell 1980). In clusters, any 
hot diffuse material originally trapped in the potentials of 
the galactic halos becomes part of the overall intraclus- 
ter medium (ICM). A galaxy whose hot gas reservoir is 
removed slowly exhausts its cold gas in a few gigayears 
(e.g. Gallagher, Bushouse, & Hunter 1989), because there 
is no supply of the fresh gas from the surrounding hot 
gas. This process is an important ingredient in the suc- 
cess of semi-analytic (SA) models, which have been shown 
to reproduce observed global properties of cluster galax- 
ies such as the morphological composition and the blue 
galaxy fraction (e.g. Baugh, Cole, & Frenk 1996; Kauff- 
mann 1996). Although such slow truncation after the in- 



fall into larger virialized halos explains the observed SFR 
gradients in clusters (Balogh, Navarro, & Morris 2000) and 
the color evolution of cluster galaxies (Kodama & Bower 
2001), the lack of the population of intermediate B/T can- 
not be solved even by including this mechanism into the 
galaxy formation model (Diaferio et al. 2001). 

The other mechanism is ram pressure stripping (RPS) of 
the cold gas from galactic disks (Gunn & Gott 1972). As a 
galaxy orbits through a cluster, it feels ram pressure (RP) 
from the ICM. When the RP is greater than the binding 
force, the cold gas will be stripped in ~ 10^ yr (Abadi, 
Moore, & Bower 1999). Fujita & Nagashima (1999) sug- 
gested that the RPS increases B/T of the Milky Way-like 
galaxies owing to the rapid suppression of the star forma- 
tion in the disk component, and consequently Sb galax- 
ies change their morphology into SO galaxies by assuming 
static cluster potentials and radially infalling galaxies. 

In this paper we investigate how these mechanisms affect 
colors and morphologies of cluster galaxies with a special 
interest in the effects of the RPS and the minor mergers 
that have not been well studied in cosmological context. 
Here we do not attempt to model the galaxy harassment, 
because we only treat bright galaxies (L > L^) and the 
galaxy harassment probably do not affect strongly for such 
bright galaxies as far as we classify the galaxies by B/T 
as we mentioned previously. 

For this purpose, the fully numerical simulations are 
still too expensive and have a lot of difficulty in handling 
subgrid physics, e.g. gas cooling, star formation, super- 
nova feedback, chemical evolution and so on, though re- 
cent work is beginning to achieve some notable successes 
for the distribution of galaxy populations (e.g. Steinmetz 
& Miiller 1995; Weinberg, Hernquist, & Katz 1997; Blan- 
ton et al. 1999; Katz, Hernquist, & Weinberg 1999; Pearce 
et al. 1999; Cen & Ostriker 2000). 

Alternative method to study the evolution of galaxy 
population is the SA modeling of galaxy formation (e.g. 
Kauffmann, White, & Guiderdoni 1993: Cole et al. 1994; 
Somerville & Primack 1999; Cole et al. 2000; Nagashima 
et al. 2001, 2002), which follows the collapse and merg- 
ing of dark halos by using a probabilistic method on the 
mass distribution based upon an extension of the Press- 
Schechter formalism (Press & Schechter 1972; Bond et al. 
1991; Bower 1991; Lacey & Cole 1993) and in which the 
subgrid physics is incorporated with the merging histories 
of dark matter halos assuming simple scaling laws. 

As a natural expansion of this approach, hybrid methods 
of cosmological A''-body simulations and the phenomeno- 
logical model of galaxy formation used in the SA tech- 
niques have been developed (Roukema et al. 1997; Kauff- 
mann et al. 1999; Benson et al. 2000). In order to ob- 
tain the spatial and velocity distribution of galaxies and 
to treat the mergers between galaxies in cluster environ- 
ment, several authors have adapted the refinements of hy- 
brid methods, which traces the merging histories of sub- 
structures within virialized halos (Okamoto & Nagashima 
2001; Springel et al. 2001). Since this method enables us to 
calculate orbits of galaxies in clusters, we can incorporate 
the RPS in our model as well as merging of substructure 
halos. 

We compare our simulated cluster at ^ = 0.2 to an ob- 
served sample that is constructed by superposing 7 low 



Environmental Effects on Cluster Galaxies 



3 



redshift CNOCl (Yee, EUingson, & Carlberg 1996) clusters 
(0.18 < z < 0.3). To compare our results with the obser- 
vations, we rescale our simulated cluster and the CNOCl 
clusters by i?2oo which is the radius of the sphere centered 
on the cluster center and whose density is 200 times the 
critical density of the universe at a given redshift. 

In §2, we describe an outline of the galaxy formation 
model used here and we show the characteristics of our 
model in §3. Results on the luminosity functions, star for- 
mation gradients, color gradients, and morphological gra- 
dients are presented in §4. These results are discussed in 
§5. 

2. MODEL 

We examine the evolution of cluster galaxies in a ACDM 

universe (JIq = 0.3, Aq = Ao/(3iJo^) = 0.7, h = Hq/WO 
km s^^ Mpc^^ = 0.7, as = 1.0). The baryon density, fib, 
is set to 0.015/1-2 (Suzuki, Yoshii, & Beers 2000). Note 
that a recent measurement of the anisotropy of the cos- 
mic microwave background suggests a slight higher value, 
^ 0.02/1-2 (Ncttcrficld ct al. 2002). The changing of 
rib affects only the normalization of parameters in our SA 
model and wc confirmed that it docs not affect our conclu- 
sion. Thus, we here prefer using fib = O.OlS/j-^ that was 
adopted in our previous paper (Okamoto & Nagashima 
2001). 

The outline of the procedures of galaxy formation is as 
follows. At first, the merging paths of dark halos are real- 
ized by a cosmological A''-body simulation. At this stage, 
it is determined for each halo whether it is a central halo 
of its virialized host halo or it is a substructure of its host. 
Next, in each merging path, evolution of the baryonic com- 
ponents, namely, gas cooling, star formation, supernova 
feedback, and chemical evolution, arc calculated based on 
the SA galaxy formation model. The detail of our SA 
model will be given in §2.3. 

We refer to a system consisting of the stars and cold 
gas as a galaxy. When two or more galactic halos merge 
together, we estimate merging time-scale of the galaxies in 
the new common halo based on dynamical friction time- 
scale. As merging of galaxies occurs, we change the mor- 
phology of the merger remnant according to the type of the 
merger. Finally, we calculate the luminosity and color of 
each galaxy based on the stellar populations that compose 
the galaxy. 

2.1. Simulations 

We first describe two simulations used in this paper and 
their specific purpose. These simulations are the same as 
those by Okamoto & Habe (1999, 2000) except for the 
background cosmology. 

For a main simulation, we adopt the constrained ran- 
dom field method to generate the initial density field in 
which a rich cluster is formed at the center of a simulation 
sphere of radius 22.5 Mpc (Hoffman & Ribak 1991). 
The constraint that we impose is the 3 a peak with the 8 
Mpc Gaussian smoothed density field at the center of the 
simulation sphere. We refer this simulation as the cluster 
simulation. 

Another simulation represents an average piece of the 
universe corresponding to the field environment within a 
sphere of radius 20 Mpc without any constraints. We 



use this simulation to normalize the parameters in our 
galaxy formation models and to compare the galaxy popu- 
lation to that in the cluster simulation. Although the size 
of this simulation seems to be somewhat small, Okamoto 
& Habe (1999) obtained a reasonable halo mass fimction 
from a simulation sphere of radius 3.5h~^ Mpc. We have 
confirmed that this simulation produces a reasonable mass 
function and the mass of the most massive halo at z = is 
sufficiently small (1.9 x IO^'^Mq) as a sample of low-density 
environment. 

To get sufficient resolution with relatively small number 
of particles we use a multimass initial condition for each 
simulation. In the field simulation, we put the high reso- 
lution particles in the sphere of radius 10 h"^ Mpc and in 
the remaining region wc put the low resolution particles as 
boundary particles. On the other hand, we adopt the res- 
imulation technique for the cluster simulation (Okamoto 
& Habe 1999). That is, first, only long- wavelength com- 
ponents are used for the realization of the initial perturba- 
tions in the simulation sphere using ~ 10^ particles, and 
then we perform a simulation with these low resolution 
particles. After this procedure, we tag the particles that 
are inside the sphere of radius 3 Mpc centered on the 
cluster center at 2 = 0. Next we go back to the initial red- 
shift, and then we divide the tagged particles according to 
the density perturbations that is produced by additional 
shorter wavelength components. As a result, the number 
of the high resolution particles becomes ^ 10^. We adopt 
the same mass and softening length of a high resolution 
particle in each simulation. Our analyses are operated 
only for these high resolution particles. 

The overall parameters and mass of the most massive 
virialized object in each simulation at 2; = are listed in 
Table 1. 

2.2. Halo Identification and Construction of Merger 

Trees 

The method to create merging history trees of dark ha- 
los is basically the same as Okamoto & Habe (2000), while 
we make some improvement. 

To construct the merger trees we identify virialized ha- 
los and their substructures at 34 redshifts during z = 20 
and z = 0. This procedure is as follows. 

At each time-step, we calculate density field by smooth- 
ing each particle with neighboring 64 particles using a 
similar method to the smoothed particle hydrodynamics 
(SPH). At a given redshift, particles with local densities 
greater than one-third of the virial density at that redshift 
are tagged. 

We then operate the friends-of-friends (FOF) grouping 
algorithm (Davis et al. 1985) only on the tagged parti- 
cles, otherwise the FOF is likely to find strangely-shaped 
halos like strings in low density regions (Bertschinger & 
Gelb 1991). The linking length is defined as bl, where I 
is the mean particle separation, and the parameter 6 is a 
function of redshift. Its value is given by 



b = 



1 Pvii 

2 p 



(1) 



where Pvir is the virial density calculated by the spherical 
collapse model at a given redshift and p is the background 
density of the universe. Note that Eq. (1) implies density 
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of a halo edge is half of the virial density that well exeeeds 
the eriterion used to tag the particles. This equation gives 
b ~ 0.2 when we assume 200 as the virial overdensity, 
which is often used in the case of the critical density uni- 
verse. When the FOF links more than 10 tagged particles, 
we identify this linked particle group as a virialized object. 
We call the halos found by this procedure FOF-halos. 

Next, we divide these FOF-halos into substructures ac- 
cording to their density peaks. For this purpose, we use 
SKID algorithm (Governato ct al. 1997) as a substruc- 
ture finder. The particles contained in the FOF-halos 
are moved along the density gradients to the local den- 
sity maxima. At each density peak, the localized parti- 
cles are grouped. After we remove the gravitationally un- 
bound particles in each group, we identify the groups that 
contain more than 10 particles as substructures. We call 
these groups SKID-halos. The SKID-halo that contains 
the most bound particle of its host FOF-halo is defined 
as a central SKID-halo and its properties are replaced by 
those of the host FOF-halo. In Fig. 1, we show the density 
map of the cluster at 2: = and the particles contained in 
the FOF-halos and the SKID-halos. 

The merger tree of each SKID-halo is constructed by the 
method described by Okamoto & Habe (2000). Note that 
we trace halo-stripped galaxies as well as the SKID-halos 
using the most bound particle in a SKID-halo as a tracer 
particle (see Okamoto & Habe 1999, 2000). 

Since cooling is not allowed in substructure halos (i.e. 
non-central SKID-halos) in our model as we will describe 
in §2.3.2, the non-central SKID-halos that do not have 
any progenitors at previous time-step cannot contain ei- 
ther stellar or cold gas components. Thus, we remove these 
useless SKID-halos from our merger trees. As a result, a 
newly-formed SKID-halo, i.e. a SKID-halo that has no 
progenitor at previous time-step, is always born as a cen- 
tral SKID-halo. 

2.3. Galaxy Formation Model 

The following prescriptions are almost the same as the 
model in Okamoto & Nagashima (2001). 

For simplicity, a dark halo is modeled as a truncated 
isothermal sphere whose mass and circular velocity are 
taken from the A''-body data. If a dark matter halo has 
no progenitor halos, the mass fraction of the gas is given 
by fh = f^b/f^o — 0.1. When a dark matter halo col- 
lapses, the gas in the halo is shock-heated to the virial 
temperature of the halo. We refer to this heated gas as 
hot gas. The gas in dense region of the halo is cooled by 
efficient radiative cooling. We dub this cooled gas cold 
gas. We estimate the amount of the gas cooled during 
a time-step using the metallicity-dcpendent cooling func- 
tion by Sutherland & Dopita (1993). In order to avoid the 
formation of unphysically large galaxies, the above cool- 
ing process is applied only to halos with circular velocities 
smaller than T4top- When we adopt the same value of 
Kitop for the field and cluster simulations, our simulations 
produce too many /few bright galaxies in the field/cluster 
simulation. Although this discrepancy is interesting, our 
simulations are too small to argue the bright end of lumi- 
nosity functions and it is beyond the scope of this paper. 
Moreover the physical meaning of l^top is quite unclear. 
Therefore we choose the values of V^top to reproduce the 



observed number of bright galaxies {L > L*). We adopt 
V^siop = 350 km s~^ and 500 km s~^ for the filed and cluster 
simulations, respectively. 

Stars are formed from the cold gas at a rate of = 
-Mcom/t*, where Mcoid is the mass of the cold gas and 
is the star formation time-scale in the disk. It means that 
the cold gas in the disk is consumed in about r, if there is 
not any supply of the fresh gas. 

When the star formation time-scale is a function of red- 
shift as in the case of t* oc Tdyn, where Tdyn is the dy- 
namical time of a galaxy, the observed properties of the 
high redshift galaxies cannot be reproduced (Kauffmann 
& Hachnelt 2000; Somerville, Primack, & Faber 2001; Na- 
gashima et al. 2001). We then adopt the description by 
Cole et al. (1994), 

where and a* are free parameters. These free param- 
eters are fixed by matching the observed mass fraction of 
the cold gas in neutral form in the disks of spiral galaxies. 

In our model, stars with masses larger than 10 Mq ex- 
plode as Type II supernovae (SNe) and heat up the sur- 
rounding cold gas. This SN feedback reheats the cold gas 
to the hot phase at a rate of Mrohcat = Afcoid/Treheat, where 
the time-scale of reheating is given by 

Treheat = (T4/T4ot)"''°*n. (3) 

where Vhot and ahot are free parameters. Although they 
are to be determined by matching the local luminosity 
function of galaxies, we set Qjiot to 2.0 because changing 
this value has a minimal effect on the results if we deter- 
mine other parameters to reproduce observations of local 
galaxies. 

Chemical evolution is treated in almost the same way 
as described in Kauffmann & Chariot (1998). The instan- 
taneous recycling approximation is adopted. We assume 
that the initially gas has a primordial composition (zero 
metallicity) . Then the gas is polluted by Typell SNe. The 
amount of metals ejected from SNe is characterized by y, 
which is heavy element yield for each generation of stars. 
The value is defined by matching the color-magnitude re- 
lation (CMR) of the local cluster ellipticals and we adopt 
y — 3Zq in this paper. All the ejecta is incorporated with 
the cold gas, while 30 % of the ejecta was assumed to be 
directly expelled to the hot gas in Kauffmann & Chariot 
(1998). The gas fraction returned by evolved stars is 0.25 
in this paper. Simultaneously, the SNe heat up the sur- 
rounding cold gas, and then metals contained in the cold 
gas are also returned to the hot gas. The returned metals 
affects cooling rate of the hot gas and the cooling transfers 
the metals from the hot gas to the cold gas. 

Before we describe our merging procedure, let us ex- 
plain our galaxy hierarchy. As we mentioned before, FOF- 
halos represent virialized dark halos. Each FOF-halo hosts 
one central SKID-halo along with a number of non-central 
SKID-halos that represent substructures of the FOF-halo. 
Each SKID-halo hosts one central galaxy along with a 
number of satellite galaxies. As we will dscribe § 2.3.2, 
only the central galaxy in the central SKID-halo can be 
supplied with the fresh cold gas by radiative cooling. The 
merging between SKID-halos are directly obtained from 
A^-body simulations. We here consider the merging be- 
tween galaxies within a SKID-halo. 
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When two or more SKID-halos have merged, we iden- 
tify a central galaxy of the largest progenitor as a central 
galaxy of the new common SKID-halo. Other galaxies are 
regarded as satellite galaxies. For simplicity, we only take 
account of the merging between a satellite galaxy and the 
central galaxy. A satellite galaxy falls into the central 
galaxy in dynamical friction time-scale, 

^ _ 1 /(^) Kr-Lo ... 
7-mrg- 2 ^ GM.atlnAc' ^ ' 

where rhaio and Vc are the radius and circular velocity of 
the new common halo, respectively, Mgat is the mass of 
the satellite, and In Ac is the Coulomb logarithm. If the 
satellite has a halo at previous time-step, Mgat is defined 
as the mass of the SKID-halo to which the satellite belongs 
as the central galaxy at the previous time-step. Since we 
trace halo-stripped galaxies as well, the mass composed by 
stars and cold gas is used as Mgat when the satellite does 
not have a halo at previous time-step. The function /(e) 
describes the dependence of the orbital decay on the ec- 
centricity of the orbit, expressed in terms of e = J/Jc{E), 
where Jc{E) is the angular momentum of a circular or- 
bit with the same energy as the satellite. This function is 
well approximated by /(e) ~ e" ''*, for e > 0.02 (Lacey & 
Cole 1993). C is a constant with value ~ 0.43. For /(e), 
we adopt the average value, (/(e)) ~ 0.5, computed by 
Tormen (1997). 

When a satellite galaxy merges with a central galaxy 
and the mass ratio of the satellite to the central galaxy 
is larger than /major, we regard this merger as a major 
merger, and then all stars of the satellite and disk stars of 
the central galaxy are incorporated with the bulge of the 
central galaxy. At the same time, cold gas contained in 
both galaxies is consumed to form bulge stars by starburst 
with the same feedback law as in the disk star formation. 

If the mass ratio does not exceed the value of the pa- 
rameter, /major, the merger is classified as minor. When 
we adopt the standard prescription, in the case of the mi- 
nor merger, the disk of the central galaxy is preserved, 
no starburst takes place, and the stars from the satellites 
are added to the bulge of the central galaxy. In this pa- 
per we also consider the case that a minor merger triggers 
small starburst because numerical experiments of merging 
galaxies have suggested that starburst is induced by a mi- 
nor merger as well as a major merger (Walker, Mihos, & 
Hernquist 1996). We dub this process minor burst and 
will give the detailed prescription in §2.3.1. 

In this paper, we adopt /major = 0.2 that was used for 
the low density universe in Okamoto & Nagashima (2001). 
We confirmed that our results hardly depend on /major, 
since our morphological classification based on B /T is nor- 
malized to reproduce the observed morphological compo- 
sition in the field. 

The absolute luminosities and colors of individual galax- 
ies are calculated using a population synthesis model by 
Kodama & Arimoto (1997). The initial stellar mass fimc- 
tion (IMF) that we adopt is the power-law IMF of Salpeter 
form with lower and upper mass limits of 0.1 and 60 M©, 
respectively. 

As we stated above, wc divided the stellar component 
into disk and bulge components. In the previous SA stud- 
ies, morphology of each galaxies have been determined by 
the B-band B/T. Simien & de Vaucouleurs (1986) showed 



that the Hubble types T of galaxies correlated with B/T 
by following relation, 

< AM >=0.324(T-^5)-0.054(T-|-5)^-|-0.0047(T-^5)^ 

(5) 

where AM = M^"^^*" - M^°**'. In this paper we use this 

relation to assign the Hubble type to our model galaxies. 
Unfortunately this relation has significant scatter (Baugh, 
Cole, & Frenk 1996), therefore wc prefer direct comparison 
of B/Ts of the model galaxies with those of the observed 
galaxies as we will describe in §4.4. 

2.3.1. Starburst triggered by a minor merger 

The effect of minor mergers was studied in Mihos & 

Hernquist (1994) using SPH simulations including star for- 
mation. In the case of a merger between a Milky Way-sized 
disk galaxy and a satellite with total mass equal to 10% of 
the central disk galaxy mass, the bar instability leads a gas 
fueling to the galactic center, and then 50 % of the total 
gas is consumed in the starburst. If the central galaxy has 
a bulge, the bulge stabilizes the disk and prevents the gas 
fueling. When the bulge-to-disk mass ratio is larger than 
1/3, the starburst becomes negligible. 

Somerville, Primack, fc Faber (2001) modeled this pro- 
cess by assuming that the fraction of cold gas consumed 
by the starburst is a function of the merging mass ratio. 
However, since only limited parameter range has been ex- 
plored, how the starburst efficiency depends on the mass 
ratio is highly unknown. Moreover, the starburst during a 
violent merger and the starburst triggered by the bar in- 
stability should be treated separate way. Then, we model 
the minor burst by assuming that the fraction, /burst, of 
the cold gas is consumed by the minor burst when the 
merger mass ratio is between /minor and /major for sim- 
plicity. In the case of the minor burst, the stellar disk 
of the central galaxy is preserved. When the bulge-to- 
disk mass ratio of the central galaxy is larger than 1/3, 
the minor burst does not take place according to Mihos 
& Hernquist (1994). We set these two parameters to re- 
produce the morphology-radius relation (/burst = 0.25 and 

/minor = 0.1/major), aS WC will show in §4.4. 

Note that if we adopt /burst = 1, it leads very similar 
results to the model without the minor burst and with 
small /major (— /minor) iu clustcrs, bccausc most of gas is 
consumed by the mergers at high redshifts in any case. 
The model with small /major only increases the bulge- 
dominated fraction and does not increase the intermediate 
B/T fraction (Okamoto & Nagashima 2001). 

2.3.2. Stripping of Diffuse Hot Gas 

We include the stripping of diffuse hot gas from galactic 
halos in our modeling as follows. When a FOF-halo con- 
tains more than two SKID-halos, we identify the central 
SKID-halo and properties of this central SKID-halo are 
replaced with those of the host FOF-halo as we mentioned 
in §2.2. The hot gas contained in the non-central SKID- 
halos are simply added to the hot gas reservoir of the host 
FOF-halo (i.e. the central SKID-halo), and so, only the 
central galaxy in the central SKID-halo can be supplied 
with the fresh cold gas by radiative cooling. This process 
has already been included in the previous SA studies and 
led significant success (e.g. Baugh, Cole, & Frenk 1996; 
Kauffmann & Chariot 1998; Diaferio et al. 2001). 
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2.3.3. Ram pressure Stripping of Cold Gas Disk 

The RPS of cold disk gas is also a new ingredient of our 
model. By following simple method, we estimate stripping 
effects of cold gas by RP from an ICM. 

In order to estimate the strength of the RP, we need 
to know orbits of cluster galaxies. For this purpose, we 
assume that the cluster has a NFW-type of density dis- 
tribution (Navarro, Frenk, & White 1997), and then we 
fit the density profile of the main cluster (i.e. the most 
massive virialized object) at each time-step of the merg- 
ing history by the NEW profile, 

'^'^ = ir/r.)illr/r.r 
where po a-nd r^ are fitting parameters. Once we fit the 
density profile of a cluster, we can calculate the poten- 
tial energy as a function of the distance from the cluster 
center, 

*(r) = -inr^Gpo^^^^^ - \HlX,r\ (7) 

Then we calculate the orbits of the galaxies within 2i?vir 
from the center of the main cluster using their positions 
and relative velocities to the cluster. 

Although gradual removal is preferable, we here consider 
only instantaneous stripping as a maximal effect because 
we have to output data with a quite small time-step to 
achieve the gradual removal of gas and such a treatment 
is not practical. Thus, the RP is estimated at pericentric 
positions for outgoing galaxies and at present positions for 
ingoing galaxies. This treatment avoids too early stripping 
for the ingoing galaxies that do not pass the pericentric 
positions during the time-step. 

The RP from the ICM is estimated as follows at the 
stripping positions, 

-Pram = PlCM^'gali (8) 

where picM is the density of the ICM and Vgai is the veloc- 
ity of a galaxy relative to the ICM. We suppose that the 
ICM distributes parallel to the dark matter for simplicity. 
The density of the ICM is obtained as j\,p(r)i where p{r) 
is taken from the NFW fit. The hot gas fraction in the 
main halo at z < 1 is 0.75, however using 0.75/bp(r) 
as picM makes no visible difference. Thus, we assume the 
hot gas fraction is unity for simplicity. 

We also attempted the ICM distribution that has the 
Coma-like core with the core density, pc = 7.98 x 10"^'' /12 
g cm~^ (Briel, Henry, & Bohringer 1992) by putting the 
ceiling to the density distribution. With the above choice 
of the core density, the core radius becomes — 140/i~^ 
kpc (~ 0.12i?2oo) &t z = 0.2. This radius is much smaller 
than the radius at which the effect of the RPS becomes 
significant as we will show in §4.2, and then we cannot 
find any difference in the results from the coreless distri- 
bution. Therefore we show only the results in the case of 
the ICM distribution parallelized to dark matter. 

For the cold gas, the restoring force per unit area due 
to the gravity of the galactic disk is given by 

-fgrav.cold = 27rGS*^disk Scold) (9) 

where S^^disk is surface density of the stellar disk and Scold 
is that of the cold gas (Gunn & Gott 1972). We cal- 
culate S*^disk as 0.5MH,_disk/7rrg, where re is an effective 



radius of the disk assuming that the effective radius is 
identical to the half-mass radius of the disk. The effec- 
tive radius is given by matching the observed luminosity- 
radius relation of spiral galaxies (see §3.1). Assuming 
Scold oc disk, the surface density of the cold gas is given 
by ■^cold — 

(Mcoid/M*,disk)S, 

.disk* 

Since we consider the maximal effect of the RPS, we 
suppose that all cold gas in a disk is stripped when 
-Pram > Pgrav • The Stripped cold gas is mixed with the 
hot gas of the host FOF-halo. 

3. MODEL SETS 

In this paper we investigate the following four baseline 
models and their combinations to clarify effect of each pro- 
cess. 

As the standard prescription of SA, we consider the 
model without either the RPS or the minor burst. By 
this model we define a reference parameter so that this 
model reproduces the properties of the local galaxies. The 
model with the reference parameter set is called the stan- 
dard model. When no fresh cold gas is supplied, the cold 
gas in a galaxy is exhausted by the star formation and 
then the star formation is truncated. Since this truncation 
time-scale is determined by the star formation time-scale, 
T*, we also consider the model with longer star forma- 
tion time-scale than the standard model. The model with 
the star formation time-scale that is four times as long as 
the standard one is referred as the low SFR model. The 
model with the minor burst and the model with the RPS 
are called the minor burst model and the RPS model, re- 
spectively. In these models, the reference parameter set is 
used. 

All models are studied in both the cluster and field sim- 
ulations except for the RPS model that is unimportant in 
the field. Thus, the corresponding model to the cluster 
standard and the RPS models in the field is the field stan- 
dard model. When we compare the cluster simulation to 
the field simulation, we always use the corresponding field 
model, for example the minor burst model in the field for 
the minor burst model in the cluster and the low SFR 
model in the field for the low SFR + RPS model in the 
cluster. These models are simply dubbed the field models. 

3.1. Param,eter setting 

In this subsection we show characteristics of our models 
and how each process affects the properties of the local 
field galaxies. The values of the parameters used in the 
standard model are listed in Table 2. 

In the top panel of Fig. 2 we show the model i?-band 
luminosity functions in the field at z = and the observed 
one taken from APM (Loveday et al. 1992), ESP (Zucca 
et al. 1997), UKST (Ratchffe et al. 1998), and 2dF (Folkes 
et al. 1999). Since the simulation volume for the field is 
too small, the brightest bin in the standard model con- 
tains only two galaxies. Therefore we do not attempt to 
reproduce the bright end of the field luminosity function 
in this paper. Apart from the bright end, our luminosity 
functions are smaller at the intermediate luminosity range 
than those in the observations and their slopes at faint- 
ends are steeper than observed slopes. This difference may 
cancel out if random collisions between satellite galaxies 
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(in substructure halos) arc considered (Somcrvillc & Pri- 
mack 1999) or if we use higher resolution that enables us 
to identify dark halos of smaller satellite galaxies (Springel 
et al. 2001). 

It is found that the minor burst hardly affects the galaxy 
luminosity function at least in the field. The luminosity 
functions in the standard and minor burst models almost 
overlap. The effect of the star formation time-scale is also 
small as shown in Cole et al. (2000). 

Although a model luminosity function is very sensitive 
to the choice of the feedback parameters (Vhot,ckhot) and 
the circular velocity above which the cooling is stopped 
(Vstop), we note that the shape of the luminosity function 
does not have particular importance for our results. The 
key observational constraints in this paper are the number 
of bright galaxies {L > L^), the cold gas mass fraction in a 
late type galaxy, and the disk size. The number of bright 
galaxies are controlled by Kstop , and latter two are used to 
specify our model parameters as described below. 

When examining the effects of gas stripping, one of the 
most important quantities is the cold gas fraction in galax- 
ies that reside in the environment where the stripping ef- 
fects arc negligible. The cold gas content of galaxies are 
strongly affected by the star formation parameters, r° and 
a* (Cole et al. 2000). In the middle panel of Fig. 2 we show 
the cold hydrogen fraction of the spiral and irregulars in 
the field models. The observational data are taken from 
the H I observation by Huchtnicicr & Richtcr (1988), which 
come from a complete sample of the galaxies with morpho- 
logical type T > 1. We classify our galaxy according to 
Eq. (5) and choose the galaxy with T > 1. It is found 
that the models with the reference parameter set repro- 
duce the observation well. There is no difference between 
the standard model and the minor burst model again. The 
fraction in the low SFR model is much larger than that in 
the standard model and the observation as expected. 

Another important quantity when we argue the RPS 
is the size of a disk that is used to calculate the surface 
density of the disk. We assume that the disk size is propor- 
tional to the virial radius of its host halo. In the bottom 
panel of Fig. 2 we show the case that re is 5 % of the virial 
radius. We also show the observed best-fit relation given 
by Totani & Yoshii (2000) for the data in Impey et al. 
(1996). It is found that our model gives a good agreement 
with the observed relation. 

One important conclusion in this subsection is that the 
minor burst hardly affects the properties of the field galax- 
ies with our values of /minor and /burst at least at ^; = 0. 

4. RESULTS 
4.1. Global Properties 

Before we investigate the distribution of galaxy popula- 
tion, it is worthwhile looking over global properties of our 
cluster galaxies. 

In Fig. 3, we show _B-band luminosity functions of clus- 
ter galaxies in the standard, low SFR, and minor burst 
models and their combination with the RPS aX z = Q. We 
also show the observed luminosity function in the Virgo 
cluster (Sandage, Binggeli, & Tammann 1985). All mod- 
els broadly agree with the observation. The luminosity 
function in the low SFR model is always larger than those 
in other models because B-band luminosity reflects recent 



star formation. The RPS affects only the low SFR model 
significantly, especially at the faint end, because the galax- 
ies in this model have more cold gas than those in the 
models with the standard star formation time-scale and 
the cold gas in the faint galaxies are easily stripped. 

To see effects of truncation of star formation on col- 
ors of galaxies, in Fig. 4 we show CMRs for all types of 
galaxies in the standard field model and galaxies within 
cluster cores (r < 0.5/i~^ Mpc) in the standard, RPS, mi- 
nor burst, low SFR and low SFR -t- RPS models. We do 
not show the minor burst -|- RPS model because it is quite 
similar to the RPS model. Here we classify the galaxies 
with T < 0.92 according to Eq. (5) as early-type galaxies 
(diamonds) and others as late-type galaxies (pluses). The 
observed CMR for the Coma ellipticals is indicated by the 
solid lines (Bower, Lucy, & Ellis 1992), which is applied 
aperture correction (Kodama et al. 1998). In the models 
with the proper star formation time-scale (standard, RPS, 
and minor burst), the blue population seen in the field 
model is vanished, and then the galaxies distribute within 
±2 mag from the CMR for the observed cluster ellipti- 
cals as observed in the Coma cluster (Terlevich, Caldwell, 
& Bower 2001). This is caused by the hot gas stripping 
in the standard model. The earlier formation epochs of 
galactic halos in the cluster environment may not be im- 
portant for this because the late-type galaxies in the low 
SFR model have significantly bluer colors than the stan- 
dard model. On the other hand, the CMRs for early-type 
galaxies are quite similar in all models, even in the field. 
It suggests that the main star formation truncation mech- 
anism for the early-type galaxies is the starburst caused 
by the major merger. 

Since the star formation is sufficiently suppressed in the 
standard model, the significant difference by the RPS is 
only seen between the low SFR model and the low SFR -|- 
RPS model. However, even in the low SFR -|- RPS model, 
there are a lot of too blue galaxies. Thus we conclude that 
the most important star formation truncation mechanism 
for the disk-dominated galaxies in cluster cores is the hot 
gas stripping. We cannot see any significant effect by the 
minor burst on the CMR. 

4.2. SFR Gradients 

Recently, many authors have analyzed radial trends in 
the SFR (Balogh et al. 1998, 1999; Elhngson et al. 2001) of 
galaxies in the CNOCl sample (Yee, EUingson, & Carlberg 
1996) . Here we investigate what mechanism is responsi- 
ble for the SFR gradient with the special interests in the 
effect of the RPS. 

In Fig. 5 we show median SFRs for observations and 
our models as a function of projected radius. The ob- 
served SFR distribution taken from Diaferio et al. (2001) 
is indicated by diamonds. The galaxies brighter than 
Mr = —20.5 are selected assuming the same cosmology 
as our simulations. In this figure we show the low redshift 
sample of Diaferio et al. (2001). The data were constructed 
by superposing 7 CNOCl clusters at the redshift range be- 
tween 0.18 < z < 0.3. Hereafter, we refer this superposing 
sample as the CNOC sample in this paper. The details 
are given in Diaferio et al. (2001) and references therein. 
We show relative SFRs, which are normalized by the field 
median value obtained from the same redshift range, be- 
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cause the observed SFRs determined by Wqu have a lot 
of uncertainties. Both for the observational and simula- 
tion data, we exclude the brightest cluster galaxies be- 
cause they often have quite different properties from other 
cluster galaxies. The standard, low SFR, and minor burst 
models are represented in the top, middle, and bottom 
panels, respectively. Corresponding RPS models are also 
shown in each panel. 

Except for the low SFR model, all our models are 
broadly consistent with the observation. As expected, the 
low SFR model has much closer SFR to the field SFR than 
other models. If we adopt the star formation time-scale 
that reproduces the observed amount of the cold gas in the 
field galaxies, the difference by the RPS is small compared 
with the observational errors. While it is interesting that 
the effect of the RPS can be seen even at r/i?2oo ~ 1-5, it 
should not be overinterpreted when comparing with obser- 
vations because the simulated region by the high resolu- 
tion particles is spherical and then the depth of the line of 
sight at outermost bin through the high resolution region 
is quite thin. As a result, in our simulation the num- 
ber of galaxies at the outer regions of the cluster becomes 
small and the field contamination is substantially under- 
estimated. For direct comparison with the observations at 
large radii, we need larger simulations. 

We confirm that the minor burst hardly affects the SFR 
distribution. 

4.3. Color Gradients 

Color of a galaxy reflects SFR of the galaxy and it is a 
more reliable observable than the SFR. It has been known 
that galaxies in inner regions of clusters are redder than 
those in outer regions (Butcher & Oemler f984). Recent 
observations have confirmed the existence of the color gra- 
dients, and besides, it continues to large radii beyond the 
virial radii (~ 2i?2oo) (Abraham et al. 1996; Carlberg et 
al. 1997). 

In Fig. 6 we show the B — V colors of the cluster galaxies 
as a function of the projected radius sk, z = 0.2. As the 
observational counterpart, we use the same sample that 
we used for the SFR gradients. Thus, we use the galaxies 
brighter than Mr = —20.5 again. The median color of 
field galaxies is shown by dotted lines. 

For the standard and minor burst models, effect of the 
RPS is more significant at large radii, because, in the cen- 
tral part, the star formation is strongly suppressed by the 
stripping of the hot gas, regardless of the RPS. Although 
it is difficult to say if the models with the RPS show bet- 
ter agreement with the observation because our models 
are less contaminated by field population as we mentioned 
previously, the standard/RPS model seems to show too 
steep/flat gradient compared with the observation. It may 
suggest that the RPS somewhat affects the color gradients 
observed in clusters. On the other hand, for the low SFR 
model, the suppression of the star formation by the strip- 
ping of the hot gas is quite weak and the cluster galaxies 
have similar colors to the field sample in this model. Then 
the suppression of the star formation is dominated by the 
RPS in the Low SFR + RPS model. It is found that the 
effect of the RPS is effective at r < 0.5i?2oo- 

While the minor burst makes galaxies slightly redder 
than the standard model, the models with and without 



the minor burst still produces very similar results. 

4.4. Morphology 

Morphologies of cluster galaxies have a strong depen- 
dence on their environments and redshifts. Here we in- 
vestigate how each process affects the morphologies of the 
model galaxies as functions of projected radius and red- 
shift. 

4.4.1. Morphology distribution 

Recently, some authors have studied the morphology- 
density or morphology-radius relations by using hybrid 
method of A'^-body simulations and SA models (Okamoto 
& Nagashima 2001; Diaferio et al. 2001; Springel et al. 
2001). Okamoto & Nagashima and Diaferio et al. con- 
cluded that the morphology evolution only by major merg- 
ers cannot explain the SO or intermediate B/T population 
in clusters and additional processes are required to under- 
stand this population. While Springel et al. seemed to 
reproduce the observational data of Whitmorc, Gilmore, 
k. Jones (1993), they adopted wider B/T range for SO's 
(0.096 < B/T < 0.45) than the standard classification 
(0.4 < B/T < 0.6). The small threshold value between 
spirals and SO's and wider range of B/T may imply neces- 
sity of the additional processes. 

One possible process is the RPS of cold gas from a galac- 
tic disk (Gunn & Gott 1972; Fujita & Nagashima 1999). 
According to Fujita & Nagashima, at a central part of a 
cluster, the SFR of the disk component rapidly drops ow- 
ing to the RPS, and then the galaxy becomes red and the 
disk becomes dark. They suggested that a Sb galaxy is 
turned into a SO galaxy as a result of the RPS. Another 
possibility is the minor merger induced bulge formation. 
Mihos & Hernquist (1994) showed that gas supply to a 
galaxy center by a minor merger triggers a star burst, and 
then the galaxy evolves toward the early-type gradually. 

We here examine how these processes affect the mor- 
phologies of cluster galaxies by comparing our galaxy pop- 
ulations in the minor burst and RPS models with those in 
the standard model and the CNOC sample. Since we do 
not have any information about galaxy morphology in our 
models except for i?/T, wc classify the galaxies based on 
their B/Ts both in our models and the CNOC sample. Ac- 
cording to Balogh et al. (1998), we classify galaxies into 
three classes: a bulge-dominated class, a disk-dominated 
class, and an intermediate class. In the CNOC sample, Di- 
aferio et al. (2001) have split the galaxy population based 
on the r-band B/T, and then galaxies with B/T > 0.7 
are identified as bulge-dominated galaxies, B/T < 0.4 as 
disk-dominated galaxies, and 0.4 < B/T < 0.7 as interme- 
diate galaxies. Wc classify our model galaxies based on the 
S-band B/T. The boundaries between the three classes 
are chosen so that the morphological fractions in the field 
models match those of the field galaxies in the CNOC sam- 
ple. As a result, the galaxies with B-band B/T > 0.4 
are identified as bulge-dominated galaxies, B/T < 0.25 as 
disk-dominated galaxies, and 0.25 < B/T < 0.4 as inter- 
mediate galaxies in the standard models, and B/T > 0.5 
for the bulge-dominated and B/T < 0.3 for the disk dom- 
inated galaxies in the minor burst models. 

The fraction of galaxies with Mn < —20.5 in each 
class is shown in Fig. 7 as a function of projected clus- 
tercentric radius and that of the field galaxies are shown 
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at r/i?200 = 2 in left panels. By comparison of the 
CNOC sample with the standard model, it is confirmed 
that the intermediate fraction in our model is much smaller 
than the observation, while the bulge-dominated fraction 
in our model agrees with the observation. This prob- 
lem has been indicated as a shortcoming of the major 
merger induced morphology evolution model (Okamoto & 
Nagashima 2001; Diaferio et al. 2001). 

When we consider the RPS, although the bulge- 
dominated fraction is increased and the disk- dominated 
fraction is decreased, the intermediate fraction is still small 
and it is even decreased at the center. It is because most 
of the disk-dominated galaxies are almost pure disk galax- 
ies as we will show and then the RPS simply darkens the 
disk-dominated galaxies rather increases their B /Ts. Fur- 
thermore the intermediate galaxies with relatively large 
B/Ts are changed into the bulge- dominated galaxies. The 
strange behavior seen in the RPS models at r/_R200 > 0.7 
is resulting from the small number of galaxies at the radius 
because of the spherical high resolution region as we dis- 
cussed previously. At that radius, all the disk-dominated 
galaxies become darker than the luminosity cut-off by the 
RPS by chance. We expect that such a behavior arc not 
observed if we consider field contamination properly and 
increase the number of clusters. In fact, when we adopt 
a fainter luminosity cut-off, this behavior is not observed 
because the number of galaxies is increased. 

The minor burst increases the intermediate fraction 
without changing the bulge-dominated fraction because it 
affects only galaxies with small bulge-to-total mass ratios. 
As a result, the minor burst model reproduces the observed 
fractions quite well in spite of it does not influence other 
properties of galaxies. Now the disk-dominated galaxies 
tend to be non-pure disk galaxies. Consequently, the RPS 
can change their B/Ts by the fading of their disks. In 
the minor burst + RPS model, we can see the increasing 
intermediate fraction towards the center. It indicates that 
the disk-dominated galaxies have larger bulges in the inner 
part of the cluster resulting from the minor burst. 

To clarify the effects of the minor burst and RPS, 
we show B/T distributions for our model galaxies with 
Mr < -20.5 at z = 0.2 in Fig. 8. The probabilistic num- 
ber distributions for the field are represented as a function 
of B/T in top panels. In the standard model, there is a 
large blank between B/T = 0.4 and 0.9, and the distri- 
bution has peaks at B/T ^ and 1. In the minor burst 
model the galaxies have similar B/T distribution to that in 
the standard model, while the peak at B/T ~ is lowered 
and the galaxies have slightly extended B/T distribution. 

In the lower panels we show the B /T distributions of the 
galaxies within -R200 iu the standard models (left) and the; 
minor burst models (right). The models with and with- 
out the RPS are indicated by the dotted and solid lines, 
respectively. To sec the influence of the RPS, we do not 
normalize the distribution by the total number of galaxies. 

In the standard model, the B/T distribution has peaks 
at B/T ~ and 1 as well as in the field model. By compar- 
ing the standard model to the RPS model, it is found that 
the RPS changes the distribution in two ways. Firstly, it 
decreases the number of galaxies brighter than the lumi- 
nosity cut-off. Since the RPS only affects the luminosity 
of the disk component, the number of the disk-dominated 



galaxies is more easily decreased than the galaxies with 
large B/Ts. Secondly, the RPS increases the B/Ts of non- 
pure disk (and of course non-pure bulge) galaxies. In our 
case, this effect is the most significant around B/T = 0.4. 

In the minor burst model, lower B/T distribution is 
significantly changed from that in the standard model. 
The peak at B/T ~ in the standard model moves to 
B/T = 0.3 — 0.4 in the minor burst model. In the cluster 
the effect of the minor burst is stronger than in the field 
because of the higher merger rate. The RPS works in the 
same way as in the standard model, but now the galaxies 
in the cluster center have larger bulges than those in the 
standard model (Fig. 7), and then the RPS increases the 
intermediate B/T fraction in the center. 

As we have shown here, the B/T distribution gives us 
much more useful information than classifying galaxies 
into three types. Therefore the direct comparison of the 
B/T distribution predicted by a model with observational 
one in various environments will be a crucial test for the 
morphology evolution model. 

4.4.2. Morphology evolution 

Dressier et al. (1997) noted that the SO fraction of clus- 
ter galaxies appears to decrease with redshift. Although 
Andreon (1998) argued that it is difficult to define a reli- 
able SO/E-ratio, van Dokkum & Pranx (2001) pointed out 
that, in clusters, at least the early-type (E+SO) fraction 
decreases with redshift. To probe the role of each environ- 
mental effect in the morphology evolution we plot the non- 
disk- dominated galaxy fractions and the number ratios of 
the intermediate galaxies to the bulge-dominated galaxies 
(hereafter Int /Bulge-ratios) in the upper and lower panels 
of Fig. 9, respectively. We select the galaxies brighter than 
Mr = —20.5 within i?vir at each redshift. 

All models show the gradually decreasing non-disk- 
dominated fractions with redshift. This qualitatively 
agrees with the observational trend of decreasing early- 
type fraction in clusters with redshift, though we cannot 
make quantitative comparison because the correlation be- 
tween the Hubble-type and B/T (Eq. (5)) has a quite 
large scatter and we do not know whether this relation is 
adequate for cluster galaxies. There is no doubt that the 
hot gas stripping and the RPS is responsible for the change 
of the non-disk-dominated fractions. These processes not 
only increase the B/T oi the galaxies but also decrease the 
number of galaxies brighter than the cut-off luminosity. 
The latter is also important because the disk-dominated 
galaxies are more easily darkened by these stripping pro- 
cesses, and then the disk-dominated fraction decreases to- 
ward lower redshift. 

Our Int/Bulge-ratios show completely different trends 
from the observed SO/E-ratio. These ratios are flat or 
moderately increasing functions of redshift, albeit the ob- 
served SO/E-ratio is a strongly decreasing function. This 
feature can be understood from the following reason. In 
our models the bulges are formed through mergers and the 
merger rate is only high before cluster formation. Thus, 
the bulges hardly grow up in clusters. Hence, the main 
process that evolves the galaxy morphology in clusters is 
the fading of disks through gas stripping. As we men- 
tioned above, this process increases B/T (i.e. it increases 
the bulge-dominated and intermediate galaxies) and dark- 
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ens the galaxies. As a result, although both the bulge- 
dominated and intermediate galaxies are increased, the 
intermediate galaxies are more likely to be darkened and 
fall below the cut-off luminosity than the bulge-dominated 
galaxies. Then the Int/Bulgc-ratios in the cluster become 
flat or increasing functions of redshift under the merger- 
driven bulge formation model. 

We have to consider the processes that form or evolve 
bulges in clusters in order to produce an Int/Bulge-ratio 
similar to the observed SO/E-ratio. However, it is proba- 
bly still possible that the Int/Bulge-ratio shows different 
evolution from the SO/E-ratio because the eyeball classifi- 
cation does not depend only on the B/T. In fact, the disks 
of luminous Sa-Sb galaxies are substantially thickened and 
the spiral features are vanished by tidal encounters in clus- 
ters and they become galaxies resembling SO's without the 
large change in their B/Ts (Moore et al. 1999). We ex- 
pect that observations of the Int/Bulge-ratio as a func- 
tion of redshift will give important information about the 
formation process and/or epoch of the intermediate B/T 
population. 

5. SUMMARY AND DISCUSSION 

In this paper, we have investigated the effects of hot 
gas stripping, RPS of cold disk gas, major mergers, and 
minor mergers on the evolution of bright cluster galaxies, 
for which the other mechanism, say, galaxy harassment, is 
probably not the predominant process. We have used a 
combination of the cosmological A'^-body simulations and 
the SA galaxy formation model. This method enables us 
to study above environmental effects in a fully cosmologi- 
cal context. We have determined the model parameters of 
the reference model to reproduce observations of nearby 
field galaxies and CMRs of local cluster ellipticals. 

Our main results are summarized as follows. 

• The process that terminates star formation in 

early- type galaxies is starburst, because the CMRs 
for early-type galaxies in the cluster are almost 
independent of either the star formation time-scale 
or the presence of the RPS. 

• If we adopt appropriate star formation time-scale, 
so as to reproduce the observed cold gas mass 
fraction in the field, the dominant process that 
determines colors of galaxies in the cluster core is 
the stripping of hot gas reservoirs. The effect of 
the RPS is only substantial when we adopt much 
longer star formation time-scale. 

• Minor burst does not affect galaxy properties 
except for morphology. 

• Without minor burst galaxies typically separate 
into almost pure disk (B/T ~ 0) and almost 
pure bulge {B/T ~ 1) galaxies. The gap in 
the B/T distribution causes the deficiency of 
the intermediate B/T galaxies and makes the 
intermediate galaxy fraction very sensitive to the 
adopted B/T range. The minor burst decreases 
the galaxies with B/T ~ and increases those 
with B/T = 0.3 - 0.4. This effect is larger in the 
cluster than in the field reflecting a richer merger 
history, and then the model with the minor burst 



can reproduce the observed morphology-radius 
relations. 

• RPS gives quite small changes in the mor- 
phological fractions in the central part of the 
cluster under the bulge formation only by major 
mergers. When the minor burst is considered, 
the disk-dominated galaxies in the cluster center 
tend to have larger bulges than in the standard 
model, so that the RPS increases/decreases the 
intermediate/disk-dominated galaxy fraction there. 

• The fractions of the non-disk-dominatcd galaxies 
in the cluster are decreasing functions of redshift. 
This feature agrees with the observed evolution 
of early-type fraction in clusters. On the other 
hand, the Int/Bulge-ratios are flat or slightly 
increasing functions of redshift unlike to the 
observed SO/E-ratio. 

The stripping of the hot gas from galactic halos sup- 
presses the star formation in the cluster core so efficiently 
that the effect of RPS is seen only at the outskirts of the 
cluster. In order to find out this effect observationally, it 
would be worthwhile comparing the colors or SFRs of the 
galaxies in cluster outskirts with those of the galaxies ly- 
ing in regions with the same local density as the cluster 
outskirts but not associated with clusters. Note that we 
do not have enough samples at the cluster outskirts and 
we do underestimate field contamination because of the 
small volume and spherical shape of the high resolution 
region in our simulation. Moreover we assume instanta- 
neous stripping of cold gas by the RP as a maximal effect, 
so that we must overestimate the stripping effect. Thus, 
it would be more diflacult to see the difference by RPS in 
observations. 

We have shown that the bulge formation solely by major 
mergers has the difficulty to form galaxies of intermediate 
B/T (~ 0.3 - 0.8) both in the field and cluster. The mi- 
nor burst can resolve this problem because it develops the 
bulge of the disk-dominated galaxies. This effect is more 
significant for the cluster galaxies than the field galaxies 
reflecting richer merger histories. Then the minor burst 
can reproduce the observed morphology-radius relation. 
Since the disk-dominated galaxies in cluster tend to have 
larger bulges than those in the field by the minor burst, 
it also gives a reasonable answer to the problem that the 
fading of the disk-dominated galaxies like field spirals can- 
not explain the entire cluster SO population as there is 
too little luminosity density in spiral bulges (e.g. Dressier 
1980). When the minor burst is considered, the RPS can 
increase/decrease the intcrmediate/disk-dominated galaxy 
fraction at the cluster center, while the colors and SFRs 
of these galaxies are hardly changed. Our results show 
that the minor burst is very important process to explain 
morphology evolution of galaxies, although our treatment 
is too simple. In order to construct a feasible morphol- 
ogy evolution model induced by mergers, the comparison 
of model predictions and observations in low density envi- 
ronments (field and groups) where the RPS and the galaxy 
harassment are not important and the merger is predomi- 
nant process that leads morphological transformation will 
be quite helpful. 
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In all our models the bulge formation is driven by merg- 
ers, and then the Int/Bulge-ratios become flat or slightly 
increasing functions of redshift contrastively to the ob- 
served SO/E-ratio in spite of we include the fading of disks 
by the hot gas and/or cold gas stripping. This is because 
the merger rate is only high before the cluster formation 
(Okamoto & Habe 1999; Gottlober, Klypin, & Kravtsov 
2001) and the bulges hardly evolve in the cluster. This 
discrepancy between our models and the observation can 
be explained by the difFcrcncc in the classification meth- 
ods, that is, the Hubble type of a galaxy depends on its 
SFR as well as its B/T. Consequently star-forming in- 
termediate B/T galaxies could be classified as spirals and 



non-star-forming disk-dominated galaxies have a chance to 
be classified as SO's. If the Int /Bulge-ratio is also observed 
as a decreasing function of redshift, we will have to con- 
sider other bulge formation processes that increase bulge 
luminosity of the disk- dominated galax;ies in high density 
environment. 

We are grateful to S. Mineshigc, Y. Fujita, and T. Ko- 
dama for their fruitful comments. Numerical computation 
in this work was carried out on SGI Origin 2000 at Di- 
vision of Physics, Graduate School of Science, Hokkaido 
University and on SGI Origin 3000 at Yukawa Institute 
Computer Facility. 



REFERENCES 



Abadi, M. G., Moore, B., & Bower, R. G. 1999, MNRAS, 308, 947 

Abraham, R. G. ct al. 1996, ApJ, 471, 694 

Andreon, S. 1998, ApJ, 501, 533 

Andrcon. S., & Ettori, S. 1999, ApJ, 516, 647 

Balogh, M. L., Scliadc, D., Morris, S. L., Yee, H. K. C, Carlberg, 

R. G., & Elliiigson, E. 1998, ApJ, 504, 75 
Balogh, M. L., Morris, S. L., Yee, H. K. C, Carlberg, R. G., & 

Ellingson, E. 1999, ApJ, 527, 54 
Balogh, M. L., Navarro, J. F., & Morris, S. L. 2000, ApJ, 540, 113 
Barnes, J. E.1996, in ASP Conference Series 92, Formation of 

the Galactic Halo — Inside and Out, ed. H. L. Morrison & A. 

Salajcdini (San Francisco: ASP), 415 
Baugh, C. M., Cole, S., & Frcnk, C. S. 1996, MNRAS, 283, 1361 
Bekki, K. 1998, ApJ, 502, L133 

Bertshinger, E., & Gclb, J. M. 1991, Comput. Phys., 5, 164 
Benson, A. J., Cole, S., Frcnk, C. S., Baugh, C. M., & Lacey, C. G. 

2000, MNRAS, 311, 793 
Blanton, M., Cen, R., Ostriker, J. P., & Strauss, M. 1999, ApJ, 522, 

590 

Bond, J. R., Cole, S., Efstathlou, G., & Kaiser, N. 1991, ApJ, 379, 
440 

Bower, R. G. 1991, MNRAS, 248, 332 

Bower, R. G., Lucy, J. R., & EUis, R. S. 1992, MNRAS, 254, 601 
Bricl, U. G., Henry, J. P., & Bohringcr, H. 1992, MNRAS, 259, L31 
Butcher, H. & Ocmlcr, A. 1978, ApJ, 219, 18 
Butcher, H. & Ocmlcr, A. 1984, ApJ, 285, 426 
Cen, R. & Ostriker, J. P. ApJ, 538, 83 

Cole, S., Aragon-Salamanca, A., Frenk, C. S., Navarro, J. F., & 

Zepf, S. E. 1994, MNRAS, 271, 781 
Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000 MNRAS, 

319, 168 

Carlberg, R. G., Yee, H. K. C, & EUingson, E. 1997, ApJ, 478, 462 
Couch, W. J., Bargcr, A. J., Small, I. Ellis R. S., & Sharpies, R. M. 

1998, ApJ, 497, 188 

Davis, M., Efstathiou, G., Frcnk, C. S., & White, S. D. M. 1985, 
ApJ, 292, 371 

Diaferio, A., Kauffmaim, G., Balogh, M. L., White, S. D. M., 

Schade, D., & Ellingson, E. 2001, MNRAS, 323, 999 
Dressier A. 1980, ApJ, 236, 351 
Dressier A. et al. 1997, ApJ, 490, 577 

EUingson, E., Lin, H., Yee, H. K. C, & Carlberg, R. G. 2001, ApJ, 

547, 609 

Folkes, S. ct al. 1999, MNRAS, 308, 459 

Fujita, Y., & Nagashima, M., 1999, ApJ, 516, 619 

Gallagher, J. S., Bushousc, H., & Hunter, D. A. 1989, AJ, 97, 700 

Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., & 

Stadel, J. 1998, MNRAS, 300, 146 
Gottlober, S., Klypin, A., & Kravtsov, A. V. 2001, ApJ, 546, 223 
Governato, F., Moore, B., Cen, R., Stadel, J., Lake, G., & Quinn 

T. 1997, New Astron., 2, 91 
Gunn, J. E., & Gott, J. R. 1972, ApJ, 176, 1 
Hoffman, Y., & Ribak, E. 1991, ApJ, 380, L5 
Huchtmeier, W. K., & Richter O.-G. 1988, A&A, 208, 237 
Impey, C. D., Sprayberry, D., Irwin, M. J., & Bothun, G. D. 1996, 

ApJS, 105, 209 

Katz, N., Hernquist, L., & Weinberg, D. H. 1999, ApJ, 523, 463 
Kauflfmann, G. 1996, MNRAS, 281, 487 
Kauffmann, G., & Chariot, S. 1998, MNRAS, 294, 705 
Kauffmann, G., Colbcrg, J. M., Diaferio, A., & White, S. D. M. 

1999, MNRAS, 397, 529 



Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 311, 576 
Kauffmann, G., White, S. D. M., & Guiderdoni, B. 1993, MNRAS, 
264, 201 

Kodama, T., & Bower, R. G. 2001, MNRAS, 321, 18 
Kodama, T., & Aiimoto, N. 1997, 

Koopmann, R. A., & Kenney, J. D. P. 1998, ApJ, 497, L751 
Kodama, T., Arimoto, N., Barger, A. J., & Aragon-Salamanca, A. 

1998 A&A 334 99. 
Lacey, C., & Cole,'s. 1993, MNRAS, 262, 627 

Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 
692 

Loveday J., Peterson, B., Efstathiou, G., & Maddox, S. 1992 ApJ, 
390, 338 

Mihos, J. C, & Hernquist, L. 1994 ApJ, 425, L13 

Moore, B., Katz, N., Lake, G., Dressier, A., & Ocmler, A., 1996, 

Nature, 379, 613 
Moore, B., Lake, G., & Katz, N. 1998, ApJ, 495, 139 
Moore, B., Lake, G., Quinn, T., & Stadel, J. 1999, MNRAS, 304, 

465 

Nagashima, M., & Gouda, N. 2001, MNRAS, 325, L13 
Nagashima, M., Totani, T., Gouda, N., & Yoshii, Y. 2001, ApJ, 
557 505 

Nagashima, M., Yoshii, Y., Totani, T., & Gouda, N. 2002, ApJ, 

578, 675 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, ApJ, 

490, 493 

Ncttcrfiel, C. B. et al. 2002 ApJ, 571, 604 
Okamoto, T., & Habe, A. 1999, ApJ, 516, 591 
Okamoto, T., & Habe, A. 2000, PAS J, 52, 457 
Okamoto, T., & Nagashima, M. 2001, ApJ, 547, 109 
Pearce et al. 1999, ApJ, 521, L99 

Press, W. H., & Schechter, P. L. 1974, ApJ, 187, 425 
RatcUffe, A., Shanks, T., Parker, Q., & Fong, R. 1998, MNRAS, 
293, 197 

Roukema, B. F., Peterson, B. A., Quinn, P. J., & Rocca- 

Volmcrrangc, B. 1997, MNRAS, 292, 835 
Sandage, A., BinggeU, B., & Tammann, G. A. 1985, AJ, 90, 1759 
Simien, F., & de Vaucouleurs, G. 1986, ApJ, 302, 564. 
Somerville, R. S., & Primack, J. R. 1999, MNRAS, 310, 1087 
Somerville, R. S., Primack, J. R., & Faber, S. M. 2001, MNRAS, 

320, 504 

Springel, V., White, S. D. M., Tormen, G., & Kauffmann, G. 2001, 

MNRAS, 328, 726 
Steinmetz, M. & Miillcr, E. 1995, MNRAS, 276, 549 
Sutherland, R., & Dopita, M. A. 1993, ApJS, 88. 253 
Suzuki, T. K., Yoshii, Y., & Beers, T. C. 2000, ApJ, 540, 99 
Terlevich, A. I., Caldwell, N., & Bower, R. G. 2001, MNRAS, 326, 

1547 

Tormen, G. 1997, MNRAS, 290, 411 

Totani, T., & Yoshii, Y. 2000, ApJ, 540, 81 

van Dokkum, P. G., & Franx, M. 2001, ApJ, 553, 90 

van Dokkum, P. G., Franx, M., Kelson, D. D., Illingworth, G. D., 

Fisher, D., & Fabricant, D. 1998, ApJ, 500, 714 
Walker, I. R., Mihos, J. C, & Hernquist, L. 1996, ApJ, 460, 121 
Weinberg, D. H., Hernquist, L., Katz, N. 1997, ApJ, 477, 8 
Whitmore, B. C, Gilmore, D. M., & Jones, C. 1993, ApJ, 407, 489 
Yee, H. K. C, Ellingson, E., & Carlberg, R. G., 1996, ApJS, 102, 

269 

Zucca E. et al. 1997, A&A, 326, 477 



Okamoto & Nagashima 



Table 1 
Parameters of Simulations 



Simulation 


Constraint 




Ni 






mh 


mi 


Rsim Mlargest 


Field 
Cluster 


None 
3(T peak 


636008 
1121534 


824438 
95406 


5e-3 
5e-3 


2e-2 
5e-2 


5.5e+8 
5.5e+8 


4.4e+9 
3.5e+10 


20 1.87e+13 ' 
22.5 5.23e+14 



Note. — Parameters A'^, e, m, Rgim^ a^nd Afiargest represent number of particles in each simulation, mass 
of a particle, radius of each simulation sphere, and mass of the largest virialized object in each simulation at 
2 = 0, respectively. Subscripts h and 1 indicate high resolution and low resolution particles, respectively. A 
length unit is Mpc and a mass unit is M© 



Table 2 

Parameters of Galaxy Formation Model 



Vhot (km s ^) 


Q!hot 


r° (Gyr) 


n, /major 


y (Zq) 


200 


2.0 


2.0 


-1.3 0.2 


3.0 
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Fig. 1. — Density map {top panel), the projection of the particles contained in FOF-halos {middle panel), and projection of those in 
SKID-halos {bottom panel) within cluster's virial radius at z = 0. 
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Fig. 2. — Properties of the galaxies in the field models as a function of i?-band luminosity. i?-band luminosity functions, cold gas mass 
fractions in spiral galaxies, and size of spiral galaxies arc presented in the top, middle, and bottom panels, respectively. The solid, dotted, 
and dashed lines indicate the standard, minor burst, and low SFR models, respectively. The observed data are represented by symbols for 
the luminosity functions (squares, asterisks, diamonds, and triangles for APM, ESP, UKST, and 2dF, respectively) and dash-doted lines for 
the cold gas mass fractions and the galaxy size. The details are given in the text. 
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Fig. 3. — -B-band luminosity functions of cluster galaxies. The diamonds, asterisks, crosses, triangles, squares, and pluses indicate the 
luminosity functions of the standard, low SFR, minor burst, RPS, low SFR + RPS, and minor burst + RPS models, respectively. The thin 
solid line is the luminosity function of the Virgo cluster galaxies by Sandage, Binggeli, & Tammann (1985). 
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Fig. 4. — Color-magnitude diagrams at 2 = 0. Wc show U — V colors of galaxies in the field model (tup left) and galaxies within the 
cluster core in the standard (top right), RPS (middle left), minor burst (middle right), low SFR (bottom left), and low SFR + RPS (bottom 
right) models as a function of V-band luminosity. Galaxies with T < 0.92 by Eq. (5) arc classified as carly-typc galaxies (diamonds) and 
other galaxies are classified as late- type galaxies (pluses). The solid lines show the observed color-magnitude relation for cluster ellipticals by 
Bower, Lucy, & Ellis (1992), which axe applied aperture correction (Kodama et al. 1998). 
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Fk;. 5. — The median SFRs of galaxies with Mjj < —20.5 arc plotted as a function of projected radius a,t z = 0.2. Solid lines in the top, 
middle, and bottom panels show the SFRs in the standard, low SFR, and minor burst models, respectively. The corresponding RPS models 
are indicated by the dotted lines. The observed median SFR of the galaxies in the CNOC sample is taken from Diaferio et al. (2001) and 
shown by the asterisks. The dot-dashed line indicates the median SFR in the field model. An error bar shows the 25th to 75th percentile of 
the distribution in each bin. 
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Fig. 6. — The median rest frame B — V colors of cluster galaxies with Mn < —20.5 are plotted as a function of projected radius. In the 
upper left panel, the diamonds denote the observed median colors of the CNOC sample, which are taken form Diafcrio ct al. (2001). The 
median color of field galaxies with Mn < —20.5 at the same rcdshift range arc shown by the dotted line. The diamonds in the upper right, 
lower left, and lower right panels indicate the median colors of the galaxies at z = 0.2 by the standard, minor burst, and low SFR models as 
a function of projected radius, respectively. The corresponding RPS model is shown by the triangles in each panel. The median color of the 
galaxies in the field model with Mn < —20.5 at 2 = 0.2 is also shown by the dotted line in each panel. An error bar shows the 25th to 75th 
percentile of the distribution in each bin. We also show the 25th to 75th percentile of the distribution for the galaxies in the field model by 
the dashed lines. 
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Fig. 7. — The morphological fractions are plotted as a function of projected radius. The diamonds, triangles, and squares indicate the 
bulge-dominated, intermediate, and disc-dominated fractions, respectively. In the top panel we show the observed fractions for the galaxies 
with Mn < —20.5 in the CNOC sample, which are taken from Diaferio et al. (2001). Other panels show the fractions in the standard, RPS, 
minor burst, and minor burst + RPS models for the simulated cluster at z = 0.2. The symbols at r/R2oo = 2 indicate the fractions in the 
field in each panel. The error bars present the Poisson errors. 



20 



Okamoto & Nagashima 



0,8 



0.6 - 



0.4 - 



0.2 - 



0,0 



— 1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



STANDARD (field) 



I ■ ■ ■ I 



0,8 



0.6 - 



0.4 - 



0.2 - 



0.0 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — 1— 



MINOR BURST (field) 



I ■ ■ . i~rn . I ... I 



0.0 0.2 0.4 0.6 0.8 
B/T 



1.0 



0.0 0.2 0.4 0.6 0.8 1.0 
B/T 




0.0 0.2 0.4 0.6 0.8 
B/T 



1.0 



20 - 



15 - 



10 - 



5 - 



1 1 1 1 — I — I — I 1 1 1 1 1 — I — I — I 1 1 1 r- 



MINOR BURST 
MINOR BURST + RPS 



I ■ ■ ■ I 



0.0 0.2 



0.4 0.6 
B/T 



0.8 1 .0 



Fig. 8. — The distribution of B-band B/T in the field (upper panels) and cluster (lower panels), we show the probabilistic distribution for 
the field galaxies and the number distribution for the cluster galaxies. In the lower panels the corresponding RPS modes are also shown by 
the dotted lines. 
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Fig. 9. — Upper. The non-disk-dominated fractions in the cluster virial radius as a function of redshift. The diamonds, triangles, squares, 
and crosses indicate the standard, RPS, minor burst, and minor burst -|- RPS models, respectively. The classification criteria axe the same as 
Fig. 7. Lower. Int/Bulge-ratio as a function of redshift. The error bars present the Poisson errors. 



